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GLOBAL TIME ESTIMATES FOR SOLUTIONS TO EQUATIONS 

OF DISSIPATIVE TYPE 

MICHAEL RUZHANSKY AND JAMES SMITH 



Abstract. Global time estimates of LP — L q norms of solutions to general strictly 
hyperbolic partial differential equations are considered. The case of special interest 
in this paper are equations exhibiting the dissipative behaviour. Results are applied 
to discuss time decay estimates for Fokker-Planck equations and for wave type 
equations with negative mass. 



1. Introduction 

The paper is devoted to the time decay of LP — L q norms of solutions to constant 
coefficients strictly hyperbolic equations of general form. It is known that such esti- 
mates lead to Strichartz estimates which are a powerful technique when dealing with 
nonlinear problems. 

We will assume that the principal part of the equation is strictly hyperbolic. The 
full equation may have variable multiplicities because of the lower order terms. One 
question of interest is to identify properties of such equations which determine the 
time decay rate of solutions. Another question of interest is what happens when there 
are multiple characteristic roots. 

Equations of higher orders appear in many applications. In particular, they arise as 
dispersion equations for hyperbolic systems, for example in the study of the Fokker- 
Planck equation and Grad systems in nonequilibrium thermodynamics. Moreover, 
in approximations of solutions to the Fokker-Planck equation the order of the cor- 
responding system tends to infinity. However, it turns out to still be possible to 
determine the decay rate of its solutions. The behaviour exhibited by these exam- 
ples is similar to the behaviour of the dissipative wave equation in the sense that 
characteristic roots lie in the complex upper half plane and come to the origin as 
single roots and at isolated points. That is why in this paper we will concentrate on 
equations of such type in Theorem I2.2[ although we will also present a more general 
Theorem 12.11 Results described here are formulated for scalar equations. However, 
they can be easily extended to systems. They also yield the well-posedness results 
for semilinear equations. Details of such analysis will appear elsewhere. 

Second order equations. 

The study of LP — L q decay estimates, or Strichartz estimates, for linear evolution 
equations began in 1970 when Robert Strichartz published two papers, [Str70aJ and 
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|Str70b| . He proved that if u = u(x,t) satisfies the Cauchy problem for the homoge- 
neous linear wave equation 



1.11 



d?u(x, t) - A x u(x, t) = 0, 0, i)eR"x (0, oo) 
u(x, 0) = 4>(x), d t u(x, 0) = ip{x), x eR n , 



where the initial data <fi and ip lie in suitable function spaces such as 
the a priori estimate 



then 



;i-2) 



{u(-,t),u t (;t),v x u{-,t))\\ Lq < c(i+ty 



||(V^,V0II^ 



holds when n > 2, p 1 + q 1 = 1, 1 < p < 2 and N p > n(p 1 



W p p stands for the standard Sobolev space with N p derivatives over L p . 
this estimate, Strichartz proved global existence and uniqueness of solutions to the 
Cauchy problem for nonlinear wave equations with suitable ("small") initial data. 
This procedure of proving an a priori estimate for a linear equation and using it, 
together with local existence of a nonlinear equation, to prove global existence and 
uniqueness for a variety of nonlinear evolution equations is now standard; a systematic 
overview, with examples including the equations of elasticity, Schrodinger equations 
and heat equations, can be found, for example, in [Rac92j . 

There are two main approaches used in order to prove (11.21) ; firstly, one may write 
the solution to (II. ip using the d'Alembert (n — 1), Poisson (n = 2) or Kirchhoff 
(n = 3) formulae, and their generalisation to large n, 



Here 
Using 



u(x, t) 



UA (2j - 1) 



dtit-% 



n — 3 

! 2 



V JdB t (x 



>dS 



n — 2 
] 2 



t" 



in-l 



Bt(x) 



(x) 
dB t (x) 



(odd n > 3) 



+(r^)^(t 



\y - x\ 



dy 



B t (x) 



•J* 



dy 



x\ 



even n) 



(for the derivation of these formulae see, for example, |Eva98] ). as is done in [vW71j 
and [Rac92j. Alternatively, one may write the solution as a sum of Fourier integral 
operators: 



u(x, t) = J 



+ e 



-it\£\ 



■0(0 



,-it\t\ 



2 ™ ' 2\£\ 

This is done in [Str70aJ, |Bre75j and |Pec76j . for example. Using one of these repre- 
sentations for the solution and techniques from either the theory of Fourier integral 
operators ( |Pec76j ). Bessel functions ( |Str70a| ) or standard analysis ( |vW71j ). the 
estimate (11.21) may be obtained. 

Another problem of interest where an LP — L q decay estimate for the linear equation 
is used to prove existence and uniqueness for the related nonlinear problem is the 
Cauchy problem for the Klein-Gordon equation. Precisely, if u = u(x,t) satisfies the 
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;i.3) 



{x,t) G W 1 x (0,oo) , 



initial value problem 

Uft(x, t) — A x u(x, t) + m 2 u(x, t) = 0, 
u(x, 0) = <p(x), u t (x, 0) = if)(x), iGl", 

where 0, if) G C^°(R ?1 ), say, and m is a constant (representing a mass term), then 

(1.4) \\(u(;t),u t (;t),V x u(;t))\\L« < C(l + t)~% \\(V x <f>, lf>)\\ <P , 

where p,q,N p are as before. Comparing (11.21) to (11.41) . we see that the estimate for 
the solution to the Klein-Gordon equation decays more rapidly — there is an improve- 
ment in the exponent of the decay function of — |(i — i). The estimate is proved 
in [vW71] . |Pec76j and [Hor97] in different ways, each suggesting reasons for this 
improvement: in [vW71j . the function 

v = v(x, x n+u t) := e- imXn+1 u(x, t) , x n+1 G R , 

is defined; using ( 11 .31) . it is simple to show that v satisfies the wave equation in 
M n+1 , and thus estimate (11.21) holds for v, yielding the desired estimate for u. This 
is elegant, but cannot easily be adapted to other situations due to the importance 
of the structures of the Klein-Gordon and wave equations for this proof. In |Pec76j 
and [H6r97j, a representation of the solution via Fourier integral operators is used 
and the stationary phase method then applied in order to obtain estimate (jl.4p . 

A third problem of interest for us is the Cauchy problem for the dissipative wave 
equation, 

u tt (x, t) - A x u(x, t) + u t (x, t) = , (x, t) G lR n x (0, oo), 
u(x, 0) = (f)(x), u t (x, 0) = if)(x), x G M n , 
G C °°(R n ). In this case, 

This is proved in [Mat 76] with a view to showing well-posedness of related semilinear 
equations. Once again, this estimate (for the solution u(x, t) itself) is better than 
that for the solution to the wave equation by — 1(^ — i); there is an even greater 
improvement for higher derivatives of the solution. As before, the proof of this may 
be done via a representation of the solution using the Fourier transform: 



where if>, 



u(x, t) 



-i 



e- 1 ' 2 sinh(| v /l-4|e| 2 ) 



V^W 2 

2e"*/ 2 sinh(|v/l-4|£| 2 ) 



+ e-*/ 2 cosh(| v /l-4|e| 2 )to) 



-i 



e-*/ 2 sin (|V4|e| 2 -l) 



+ e 



-t/2 



cos(|v/4|e| 2 -l) 0(0 



2e-'/ 2 sin(|v/4|e| 2 -l) 



m), iei>i/2. 



Matsumura divides the phase space into the regions where the solution has different 
properties and then uses standard techniques from analysis. 
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Problem. 

It is, therefore, interesting to ask why the addition of lower order terms improves the 
rate of decay of the solution to the equation; furthermore, we would like to understand 
why the improvement in the decay is the same for both the addition of a mass term 
and for the addition of a dissipative term. In the proof of each of the estimates 
(see the papers cited above), the critical role is played by the characteristic roots 
of the equations. In fact, it is the difference in the behaviour of the characteristic 
roots of the Klein-Gordon equation and the dissipative wave equation which yield 
improvement over the decay rate for the wave equation. 

The aim of this paper is to investigate this phenomenon for higher order hyperbolic 
equations and see how lower order terms affect the rate of decay compared to that 
for the homogeneous m th order equation and the examples above. Equations of this 
type appear in many applications. In particular, they arise as dispersion equations 
for hyperbolic m x m systems. The order m may be large, as in, for example, Grad 
systems coming from nonequilibrium thermodynamics, where it corresponds to the 
number of moments under consideration. Moreover, in applications to Fokker-Planck 
equations describing the distribution of Brownian particles, the order m corresponds 
to the Galerkin approximation of solutions, so it is increasing to infinity. In all these 
cases equations become too large and involved to analyse explicitly, so we are led to 
study properties which determine the decay rate of L p — L q estimates in the general 
form. 

We will consider the Cauchy problem for m th order constant coefficient linear 
strictly hyperbolic equation of the general form for u = u(x,t): 

{m m—1 
D T U + E P j( D x)D™~ j u + ^2 E c a , r D a x Dlu = 0, t>0, 
3=1 1=0 \a\+r=l 

D l t u(x, 0) = ft{x) e C™(R n ), I = 0, . . . , m - 1, xeW 1 , 

where Pj(£) is a constant coefficient homogeneous polynomial of order j, and the c a , r 
are constants. 

We seek a priori estimates for the solution to this problem of the type 

m—1 

(1-6) \\D%Dlu{;t)\\L« < K{t) , 

1=0 P 

where l<p<2, - + - — 1, N p — N p (a,r) is a constant depending on p, a and r, 
and K(t) is a function to be determined. 

Homogeneous equations. 

The case where the operator in (11.51) is homogeneous has been studied extensively 
and provides many interesting relations to the geometric properties of characteristics. 
In this case we have 



(1.7) 



L(D X , D t )u = 0, (x, t) e W n x (0, oo) 
D[u{x, 0) = fi(x), I = 0, . . . , m - 1, x e W 1 , 
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where L is a homogeneous m th order constant coefficient strictly hyperbolic differen- 
tial operator; the symbol of L may be written in the form 

L(t, = (r - ^(0) ■ • ■ (r - <p m (0), with ^(0 > ■ ■ ■ > <p m (£) (f ^ 0). 

In a series of papers, |Sug94| , |Sug96| and |Sug98| , Sugimoto showed how the geomet- 
ric properties of the characteristic roots ipi(£), ■ ■ ■ , <Pm{Q affect the LP — L g estimate. 
To understand this, let us summarise the method of approach. 

Firstly, the solution can be written as the sum of Fourier multipliers: 

m— 1 m 

u(x,t) = J2lEi(t)fi](x)i where E t (t) = ^ ^V^a^O? 

1=0 k=l 

and Ofc,z(£) is homogeneous of order —I. Now, the problem of finding an LP — LP decay 
estimate for the solution is reduced to showing that operators of the form 

M r (D) ^-v^rxto?, 

where G C""(R n \ {0}) is homogeneous of order 1 and \ G C°°(W L ) is equal to 1 
for large £ and zero near the origin, are LP — LP bounded for suitably large r > I. In 
particular, this means that, for such r, 

\mi)f\\ Lq <c\\f\\ wrl . 

Indeed, it may be assumed, without loss of generality, that t — 1 since for t > 
and / G C^°(R n ), we easily have 

[E l (t)f](x)=t l [E l (l)f(t.)](r l x). 

Using this identity gives 

\\w)f\\h=*\\[W)ft)(t- i Wr*=i! q / imrnKt-'x^dx 

f r|[^(i)/ t ](^)i 9 ^' = ^ + 1l^(i)/tllL- 



w k — \\J \\w k ' 



Then, noting that a simple change of variables yields 

we have, 
hence, 



\\ft\\ P Wpk < c& 

mt)f\\ Lq <ct i+ i\\f t \\ wrl <cf- n ^ ■ 



m—l 

■ l l\ 



\\u{-ML«<c?- n ^Yj^\ wi -i. 

1=0 

It has long been known that the values of r for which M r (D) is LP — LP bounded 
depends on the geometry of the level set 

E v = {£ G R n \ {0} : <p(£) = 1} . 

In |Lit73j . |Bre75j it is shown that if the Gaussian curvature of is never zero 
then M r (D) is LP — L q bounded when r > ^^(^ — This is extended in |Bre77j . 

where it is proven that M r (D) is L p — LP bounded provided r > ^-^(^ — -), where 
p = min^ rank Hess <£>(£). 
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Sugimoto extended this further in |Sug94] , where he showed that if E<p is convex 



then M r (D) is L p - L q bounded when r > (n - ^4) (- - -); here 



7(E) := sup sup 7(E; cr,P) , EcK"a hypersurface , 
o-es p 

where P is a plane containing the normal to E at o and 7(S; a, P) denotes the order 
of the contact between the line T a R P, T a is the tangent plane at <r, and the curve 

£n P. 

In order to apply this result to the solution of (11.71) . it is necessary to find a condition 
under which the level sets of the characteristic roots are convex. The following 
notion is the one that is sufficient. Let L = L(D X , D t ) be a homogeneous m th order 
constant coefficient partial differential operator. It is said to satisfy the convexity 
condition if the Hessian, Hess<^fc(£), corresponding to each of its characteristic roots 
<fi(C), ■ ■ ■ 1 (Pm(0 is semi-definite for £ 7^ 0. 

It can be shown that if an operator L does satisfy this convexity condition, then 
the above results can be applied to the solution and thus an estimate of the form 
(Oil holds with 

Kit) = (i + ty^ip-*) , 

where 7 = maxi<fc< m 7(S^). We also have 7 < m. 

Finally, in the case when this convexity condition does not hold, it was shown 
in |Sug96| and Sug98 that, in general, M r (D) is LP — LP bounded when r > (n — 



70 (£) 



where 



70(E) := sup inf 7(E; a, P) < 7(E). 

(TEE P 



For n = 2, 70(E) = 7(E), so, the convexity condition may be lifted in that case. 
However, in |Sug96|, examples are given when n > 3, p = 1, 2 where this lower bound 
for r is the best possible and, thus, the convexity condition is necessary for the above 
estimate. It turns out that the case n>3, l<p<2is more interesting and is studied 
in greater depth in |Sug98|, where microlocal geometric properties must be looked 



at in order to obtain an optimal result. It can be noted that in L p — L p estimates 
other geometric properties of phase function and wave fronts become important, see 
the survey [RuzhOOj for more details. 

Two remarks are worth making; firstly, the convexity condition result recovers 
the Strichartz decay estimate for the wave equation, since that clearly satisfies such 
a condition, Secondly, the convexity condition is an important restriction on the 
geometry of the characteristic roots that affects the LP — L q decay rate; hence, in the 
case of an m th order operator with lower order terms we must expect some geometrical 
conditions on the characteristic roots to obtain decay. 



2. Main Results 



In this paper we will present conditions under which we can obtain LP — LP decay 
estimates for the general m th order linear, constant coefficient, strictly hyperbolic 
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Cauchy problem 

(2.1) ' ' 

{m m—l 
P(D t ,D x ) = D™u + Y,Pj(D x )D™->u + J2 Yl c a ,rD a x D r t u = Q, t > 0, 
3=1 1=0 \a\+r=l 

D l t u(x, 0) = fi{x) G C?(R n ), I = 0, . . . , m - 1, x G W 1 . 

As usual, the strict hyperbolicity means that the principal symbol of the operator 
P(D t , D x ) is strictly hyperbolic, i.e. has real roots, distinct for £ 7^ 0. However, 
characteristic roots T\ (£),..., r m (£) of the full symbol may have any multiplicities. 
Since we are interested in the question of how do lower order terms influence time 
decay rates, we do not want to worry about the well-posedness of the Cauchy problem 
and therefore assume that the principal part of the operator is strictly hyperbolic. 
Our main Theorem 12.11 states how different behaviour of the characteristic roots 
T\ (£),..., r m (£) affect the rate of decay that can be obtained. For now we will as- 
sume that symbols of Pj(D x ) are homogeneous polynomials of order j with constant 
coefficients. However, results extend to the case when Pj are pseudo-differential op- 
erators of order j and when lower order terms are pseudo-differential in D x , provided 
the statement of Lemma [3^1 holds. This case is essential when considering hyperbolic 
systems and their dispersion equations. 
It is natural to impose the condition: 

(2.2) lmr fc (0 > for k = 1, . . . , m and for all £ G W l ; 

this is equivalent to requiring the characteristic polynomial of the operator to be stable 
at all points £ G R™, and thus cannot be lifted, since we can not expect any time decay 
if this condition fails. Coefficients of equation (12. ip are allowed to be complex as long 
as condition (12. 2p holds. Of course, because of the strict hyperbolicity, coefficients of 
the principal part are real. 

Also, it is sensible to divide the considerations of how characteristic roots behave 
into two parts: their behaviour for large values of |£| and for bounded values of 
These two cases are then subdivided further; in particular the following are the key 
properties to consider: 

• multiplicities of roots (this only occurs in the case of bounded |£|); 

• whether roots lie on the real axis or are separated from it; 

• behaviour as |£| — >• 00 (only in the case of large |£|); 

• how roots meet the real axis (if they do) ; 

• properties of the Hessian of the root, Hessrfc(£); 

• a convexity-type condition, as in the case of homogeneous roots. 

Some definitions will be needed for the main theorem. Given a smooth function 
r : R n -> R and A G R, set 

S A = E A (r) := {( 6 K™ : r(£) = A} . 

In the case where r(£) is homogeneous of order one, write S T := Si(r) — for such r, 
we then have S A (r) = AS T . Also, a smooth function r : W 1 — > M will be said to 
satisfy the convexity condition if S A is convex for each A G I. Note that the empty 
set is considered to be convex. Finally, we will use notation in the introduction for 
the maximal orders of contact of a hypersurface, 7(E) and 70(E). We note that 
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if is a polynomial of order m and S = {£ G IR n : = 0} is compact then 
7o(S) < 7(S) < m; this is useful when applying the result below to hyperbolic 
differential equations and is proved in |Sug96|. 
Now we may state the main theorem: 

Theorem 2.1. Suppose u = u(x,t) satisfies the m th order linear, constant coeffi- 
cient, strictly hyperbolic Cauchy problem 112.1}) . Denote the characteristic roots of the 
operator by Ti(£), . . . , r m (£) and assume that 112.2}) holds. 

We introduce two functions, K^(t) and K^ b \t), which take values as follows: 
I. Consider the behaviour of each characteristic root, 7fc(£), in the region |£| > N , 
where N is some large number. The following table gives values for the function 
Kfr(t) corresponding to possible properties of if Tfc(£) satisfies more than 
one, then take K^(t) to be function that decays the slowest as t — > oo. 



Location of Tfc(£) 


Additional Property 




away from real axis 




e~ dt , some <5 > 


on real axis 


detHessr fc (0 ^ 
rank Hess (£) = n — 1 
convexity condition, 7 
no convexity condition, 70 


(l+t) a(„ ,) 

U — 1/1 1\ 

(1 + t) ~> 9 } 

, n-1/1 In 
(1 + t) ~ ( f 1> 
1/1 In 
(l+t) Vp *' 


asymptotic to real axis 


detHessr fc (0 ^ 
rank Hess (£) = n — 1 
no convexity condition, 70 


(l+t) 2 ( f ^ 
, n-1/1 In 
(1 +t) ~> « j 
1/1 In 

(l+t) Vp 1' 



Then take K (l) (t) = max fc= i ..., m j (t) . 

II. Consider the behaviour of the characteristic roots in the bounded region |£| < AT; 
again, take K^ h \t) to be the maximum {slowest decaying) function for which 
there are roots satisfying the conditions in the following table: 



Location of Root(s) 


Properties 


K (h \t) 


away from axis 


no multiplicities 
L roots coinciding 


e~ At , some S > 
(1 + t) L ~ l e- St 


on axis, 
no multiplicities 


detHessT fc (f) ^ 
convexity condition, 7 
no convexity condition, 70 


(l+t) 2 ( f i ] 

n I In 
(1 + t) 1 { P l' 
1,1 In 

(l + t) Vp i> 


meeting axis 
with finite order s 


L roots coincide 
on set of codimension i 


(1 + 1) SV P 1' 



Then, with K(t) = max (K^(t), K^(t)) , the following estimate holds: 

m— 1 

\\D a x Dlu{;t)\\ Lq <K(t)Y,\\fi 



1=0 



where 1 < p < 2, - + ~ = 1, and N p = N p (a, r) is a constant depending on p, a and r. 

Let us make a number of remarks on how to understand this theorem. Since the 
decay rate does depend on the behaviour of characteristic roots at different points, 
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we single out properties which determine this decay rate. Since the same character- 
istic root, say t&, may exhibit different properties at different points, we look at the 
corresponding rates K^°\t), K®(t) under each possible condition and then take the 
slowest one for the final answer. It also means that if we microlocalise in a region 
where only one of these properties holds, we can get the decay rate straight from the 
table for the corresponding solution. In some cases, especially when roots do not lie 
on the axis for large £, the result may be extended to p = 1. 

In Part I of the statement, it can be shown by the perturbation arguments that 
only three cases are possible for large £, namely, the characteristic root may be 
uniformly separated from the real axis, it may lie on the axis, or it may converge to 
the real axis at infinity. If, for example, the root lies on the axis and, in addition, it 
satisfies the convexity condition with index 7, we get the corresponding decay rate 

K®(t) = ~r^~p~Q . Indices 7 and 70 in the tables are defined as the maximum 

of the corresponding indices 7(£a) and 7(Xa), where £a = {£ : Tfc(£) = A}, over all k 
and over all A, for which £ lies in the corresponding zone. 

The statement in Part II is more involved since we may have multiple roots in- 
tersecting on rather irregular sets. The number L of coinciding roots corresponds to 
the number of roots which actually contribute to the loss of regularity. For example, 
operator (df — A) (<9 t 2 — 2 A) would have L = 2 for both pairs of roots intersecting at 
the origin. Meeting the axis with finite order s means that we have the estimate 

(2.3) dist(£,Z fc ) s <c|Imr fc (0| 

for all the intersecting roots, where Z}~ = {£ : Im Tfc(£) = 0}. In Part II of Theorem l2.ll 
the condition that L roots meet the axis with finite order s on a set of codimension £ 
means that all these estimates hold and that there is a (regular) set Z of codimension 
£ such that Z^ C Z for all corresponding k. In Theorem 12.21 we will discuss the 
special case of a single root meeting the axis at a point £0 with order s, which 
means that Imrfc(£o) = and that we have the estimate |£ — £o| s < c |I mT fc(OI- I* 1 
fact, under certain conditions an improvement in this part of the estimates is possible, 
see Theorem 12.21 and Remark 12.31 

In addition to the theorem, if we have L multiple roots which coincide on the real 
axis on a set S of codimension £, we have an estimate 

m— 1 

(2-4) \u(t,x)\ ^Cil + t^'^WMW, 

1=0 

if we cut off the Fourier transform of the Cauchy data to the e-neighbourhood S e 
of S with e = 1/t. Here we may relax the definition of the intersection above and 
say that if L roots coincide on a set S, then they coincide on a set of codimension 
£ if the measure of the e- neighborhood S e of S satisfies \S € \ < Ce e for small e > 0; 
here S' = e R" : dist(£, S) < e}. The estimate (I2.4p follows from the procedure 
described below of the resolution of multiple roots. We can then combine this with the 
remaining cases outside of this neighborhood, where it is possible to establish decay 
by different arguments. In particular, this is the case of homogeneous equations with 
roots intersecting at the origin. However, one sometimes needs to introduce special 
norms to handle L 2 -estimates around the multiplicities. Details of this will appear 
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elsewhere. Finally, in the case of a simple root we may set L — 1, and I — n, if it 
meets the axis at a point. 

Theorem 12.11 allows a microlocalisation and estimates for the corresponding oscil- 
latory integrals. In fact, Theorem 12.11 follows from its microlocal version in regions 
where characteristic roots are simple. In regions with multiple roots one requires 
additional arguments resolving the singularities caused by multiple roots followed by 
estimates for relevant pieces of the solution. A microlocal version of the theorem leads 
to better estimates since in Theorem 12.11 we take the slowest among all microlocal 
decay rates. Such a microlocal version and the full proof of Theorem 12.11 will appear 
elsewhere. 

For our applications in Section HJ we only need this result in the special case where 
characteristic roots meet the real axis with finite order; therefore, we shall state and 
outline the prove the theorem in this special case. 

Theorem 2.2. Consider the m th order strictly hyperbolic Cauchy problem A2.1\) for 

r»i_i-i4.| Q (l_j4- r 

operator P(D t , D x ), with initial data fj £ L p fl W^ 2 for j — 0, . . . , m — 1, 

where 1 < p < 2 and 2 < q < oo are such that - + - = 1. Assume that the 

— f — — H — V I 

characteristic roots 7~i(£), . . . , r m (£) of P(r, £) = satisfy \2.2\) and the following 
conditions: 

(HI) there is some e > such that for all k £ {1, . . . , m} we have 



(H2) for each £ G M n there is at most one k for which ImTfc(£o) = and there exists 
a constant c > such that 



for £ in some neighbourhood of £o- 
Then the solution u = u(x,t) to Cauchy problem \2. 1\) satisfies the estimate 



Essentially, this theorem is a special case of Theorem 12.11 where we get the ex- 
ponential decay from Part I, exponential decay from multiple roots away from the 
real axis in Part II, as well as the last line of the table in Part II with L = 1 and 
I = n, since we have only a single root coming to the axis. The main problem here is 
the possible appearance of multiple roots in the complex upper half plane. If several 
roots meet on the axis, the decay is then given in Part II of Theorem 12.11 where we 
observe the appearance of the extra power t L ~ x compared to Theorem 12.21 If roots 
come to the axis on a set of other codimension £, the order should change according 
to Theorem 12.11 If conditions of Theorem 12.21 hold only with £ we will call the poly- 
nomial P(t, £) strongly stable. Such polynomials will be discussed in more detail in 
Section HI 



lim inf Im rfc (£) > £ ; 



|£-£o| s <c|lmr fc (0|, 



m— 1 
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Remark 2.3: The order of time decay in Theorem 12.21 may be improved in the 
following cases. If Imr fe (^ ) = in (H2) implies that £ — 0, then we actually get 

m—l 

Now, assume that for all £o in (H2) we also have the estimate 
(2.5) c |e - U s < c|lmr fc (0| < c x |£ - U°\ 

with some constants Co, c\ > 0. 

If Imr/c^o) = in (H2) implies that ReT^o) = 0, then we actually get 

m— 1 



D r t D>(;t) t <C(l + t)-(^)(^)j2\\f : 



j=0 w p v VLYJ ; 



And finally, assume that for all £o such that Imrt((o) = in (H2), we also have 
£o = and Rerfc(^ ) = 0. Then we actually get 



L?(K«) 



m-1 

<c(i + ty' ±h ^ 1 ( l P -l)y m 2 , , 



+ |a|+r — j 



The proof of this is based on Remark 13.81 In particular, these estimates cover the 
case of dissipative wave equation and applications in Section HI A similar remark can 
be made for Theorem 12. 1[ where we also get the corresponding improvements. 



3. Outline of the proof 

Here we will outline the proof of Theorem 12.21 For large frequencies we have 
simple roots separated from the real axis so we can expect exponential decay in time 
there. For small frequencies, while separated from the real axis, we may have multiple 
roots, which may intersect on a rather irregular set. We will cut off around this set 
and show that we can get additional polynomial growth in time dependent on the 
"dimension" of this set, which is matched against exponential decay. Technically we 
have to establish a number of additional estimates on the solution in this case since 
the usual solution representation blows up around points of multiplicity. Finally, we 
can show the polynomial decay in time when characteristic roots approach the real 
axis. 

3.1. Some properties of hyperbolic polynomials. Here we will describe some 
useful properties of hyperbolic polynomials. Let L = L(D X , D t ) be a linear m th order 
constant coefficient partial differential operator. Then each of the characteristic roots 
of L, denoted Ti (£),... , T m (£), is continuous in R 71 ; furthermore, for each k — 1, . . . , m, 
the characteristic root Tfc(£) is analytic in 

Let now L = L(D x ,D t ) be a linear m th order constant coefficient strictly hyperbolic 
partial differential operator. Then there exists a constant iV such that, the charac- 
teristic roots Ti(£), . . . , r m (£) of L are pairwise distinct for |£| > N. We also have the 
following symbolic properties of characteristic roots: 
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Proposition 3.1. Let L = L(D x ,D t ) be a linear m th order constant coefficient hy- 
perbolic partial differential operator with characteristic roots ri(£), . . . , r m (£); then 

I. for each k = 1, . . . , m, there exists a constant C > such that 

\T k (Z)\<C(l + \Z\) for all eel n . 

Suppose that the maximum order of the lower order terms is < K < m — 1. 
Furthermore, assume that L is strictly hyperbolic, and denote the roots of the principal 
part L m (£,r) by <fii(l;), ■ ■ ■ ,y?m(0- Then we have the following: 

II. For each r k (^), k = 1, . . . ,m, there exists a corresponding root of the principal 
symbol (possibly after reordering) such that 

\T k (0-MO\<C(l + \Z\) K+1 - m for all {eW 1 . 

III. There exists N > such that, for each characteristic root of L and for each 
multi-index a, we can find constants C = C k , a > such that 

\d?r k (0\ <C\^-\ a \ forall\i\>N, 

IV. There exists N > such that, for each T k (£,) a corresponding root of the 
principal symbol can be found (possibly after reordering) which satisfies, 
for each multi-index a and k = 1, . . . , m, 

|a«r fe (0 - d«MO\ < C\i\ K+1 - m -^ for all |d > N 

for each multi-index a and k = 1, . . . , m. 

3.2. Representation of the solution. Recall that we begin with the Cauchy prob- 
lem with solution u = u(x,t): 



(3.1) 



m rn—1 

D?u + Y,PADz)DT~ 3 u + Y, E c a , r D a x D r t u = 0, t > 0, 

j=l 1=0 \a\+r=l 

D[u(x, 0) = fi(x) G L7 °°(R n ), I = 0, . . . , m - 1, x G W l , 



where symbol Pj(£) of Pj(D) is a constant coefficient homogeneous polynomial of 
order j, and the constants. 

Applying the partial Fourier transform with respect to x yields an ordinary differ- 
ential equation for u = tt(£, t) := L„ e~ tx '^u(x, t) dx: 

m m—1 

(3.2a) D?u + J2Pj(Z)Dr j u + J2 E c a ,rCD r t u = 0, 

3=1 1=0 \a\+r=l 

(3.2b) D l t u(£,0) = f l (O, l = 0,...,m-l, 

where (£, t) G W 1 x [0,oo). Let Ej = Ej(£,t), j = 0,...,m — 1, be the solutions 
to f)3.2ap with initial data 

(3.2c) D l t E^,0) = 
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Then the solution u of (13. II) can be written in the form 

m— 1 

(3.3) u(x,t) = Y,(?~ 1 E J ?f 3 ){x,t) 1 

j=0 

where 7 and represent the partial Fourier transform with respect to x and its 
inverse respectively. 

Now, as (I3.2al) . (I3.2cj) is the Cauchy problem for a linear ordinary differential 
equation, we can write, denoting the characteristic roots of ( 13. ip by Ti (£),... ,T m (£), 



E^,t) = J2 A i^ e 



»7* (0* 

"-j vsi 

fc=i 



where Aj(£,t) are polynomials in t whose coefficients depend on £. Moreover, for 
each k = 1, . . . , m and j = . . . , m — 1, the Aj(£, t) are independent of t at points 
of the (open) set {( 6 1" : T k(Q r z(0^ k}i when this is the case, we write 
Aj(£,t) = Aj(£). For Aj(£), we have the following properties: 

Lemma 3.2. Suppose ( 6 := {( G 1" : r fc(0 7^ T i{0^^ ^ en u ' e ^oi>e the 

following formula: 

m—j—l 

E fc n ^(o 

(3.4) 4(0 = ^<^-^ , 

II (n(0-r k (0) 

l=l,l^k 

where J^ fe means sum over the range indicated excluding k. Furthermore, we have, 
for each j — 0, . . . , m — 1 and k — 1, . . . , m, 

(i) Aj(£) is smooth in 5&; 
(«) Aj(£) = 0(\t\-*) as |£|-oo. 

Proof. The representation (13 ,4p follows from Cramer's rule (and is done explicitly 
in |Kli67] ): = , where V := (r^ -1 (£))"!_.. is the Vandermonde ma- 

trix and Vlf is the matrix obtained by taking V and replacing the k th column by 
( .. . 1 . . . 0) T . 

3 

Smoothness of A 1 ?^) in Sk is obvious and the asymptotic behaviour is a consequence 
of lPart II of Proposition O since (E3D holds for all |f | > N. □ 

In view of Lemma l3~2l choose N\ > so that the Tfc(£), k — 1, . . . , n, are distinct for 
|£| > iVi. Also, choose iV2 > so that all points at which any of the roots, Tfc(£), meet 
the real axis — i.e. points £ G R n such that, for all e > 0, there exist £1,^2 £ ^e(0 
with Imr fc (6) = and Imr fc (£ 2 ) 7^ 0— lie in B N2 (0). Set JV = max(JVi, JV 2 ). 

Let x(0 = Xn(0 e C^°(R n ), < x(0 < 1, be a cut-off function that is identi- 
cally 1 for |£| < iV and identically zero for |£| > 2iV. Then (13.31) can be rewritten 
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as: 



m—l 



m—l 



(3.5) u(x,t) = J2^\E, X 3 : mx,t) + J2^ 1 (E j (l- X )3 : f j )(x,t) 

3=0 3=0 



3.3. Large |£|. The second term of ( 13.51) is the most straightforward to study: by 
the choice of N, 



^(u)(i-x)(o = E4(0(i-x)(^' 



fe=l 

therefore, since each summand is smooth in W 1 , 

m—l 



Y J ^- 1 {E j {l- X )Jf j ){x 1 t) 



3=0 



m—l rn 



n > j=0 fc=i 



(2 

Note that, unlike in the case of homogeneous strictly hyperbolic equations we may 
not assume that £ = 1. Each of these integrals may be studied separately. Indeed, 
we have the following result: 

Proposition 3.3. Let r : U — > C be a smooth function, U C M n open, and dj G 
Siq(U). Assume: 

(i) there exists 5 > such that Imr(£) > 8 for all £ G U; 

(ii) \t{Z)\<C{1 + \£\) for aH£e U. 
Then, 



%(0fV(07(0^ 



< Ce 



-<St 



r .a= 1" 



If 



JV + |a|+r-j 



and 



L2(RS) 



< Ce 



-<5i 



/or a// £ > 0, iVo > n, multi-indices a, r G M and / G C^°(U). 
So, for all £ > 0, 



££££ 



3 i(a:-M-r(£)i) 



aj(0f(0 dx 



< Ce 



-St 



w N 1 + \a\+r-j , 



3 i(x-f+r(Ot) 



L- 



< Ce 



-at 



if. 



l"l + >— j • 



where iVi > n, r > 0, a multi-index; interpolating then gives, 



did: 



aj(0f(0dx 



< Ce 



-51 



Li 



where p 1 + q 1 = l, l<p<2, N p > n(- — -), r > 0, a a multi-index and 
/ G C^°(R Tl ). Thus, in this case we have exponential decay of the solution. 



GLOBAL TIME ESTIMATES FOR SOLUTIONS TO EQUATIONS OF DISSIPATIVE TYPE 15 



3.4. Bounded |£|. Let us now consider the terms of the first sum in (|3.5p . the case 
of low frequencies, 

i r m 

(3.6) ?-\E jX ?f)(x,t) = - / e^(^e^^4(e,t))x(0/(0^. 

Unlike in the case above, here the characteristic roots Ti (£),..., r m (£) are not nec- 
essarily distinct at all points in the support of the integrand (which is contained in 
the ball of radius 2N about the origin); in particular, this means that the A}(f,t) 
genuinely depend on t and we have no simple formula valid for them in the whole 
region. 

For this reason, we begin by systematically separating neighbourhoods of points 
where roots meet — referred to henceforth as multiplicities — from the rest of the re- 
gion, and then considering the two cases separately. 

First, we need to understand in what type of sets the roots 7>.(£) can intersect: 

Lemma 3.4. The complement of the set of multiplicities of a linear strictly hyperbolic 
constant coefficient partial differential operator L(D X , D t ), 

S := {£ G R n : r,(0 ^ r k (0 for all j ^ k} , 

is dense in M. n . 
Proof. First note 

S = {£eR n : A L (0 ^ 0} , 

where A^ is the discriminant of L(£, r). Now, by Sylvester's Formula (see, for exam- 
ple, [GKZ94J), Al is a polynomial in the coefficients of L(£,r), which are themselves 
polynomials in £. Hence, A^ is a polynomial in £; as it is not identically zero (for 
large |£|, the characteristic roots are distinct, and hence it is non-zero at such points), 
it cannot be zero on an open set, and hence its complement is dense in MJ 1 . □ 

Corollary 3.5. Let L(£,r) be a linear strictly hyperbolic constant coefficient partial 
differential operator with characteristic roots Ti(£), . . . , r m (£). Suppose M. k i C W 1 is 
a set such that Tfc(£) = tj(£), for some k ^ I, for all £ G M-u- For e > 0, define 

Ml l :^{Ce^ n --dist(C,M kl )<e}; 

denote the minimal v G N such that meas(A^| ; ) < Ce v for all sufficiently small e > 
by codimyWfci. Then codimA^fc/ > 1. 

Proof. Follows straight from Lemma [3.41 the fact that M. k i has non-empty interior 
ensures that its ^-neighbourhood is bounded by Ce in at least one dimension for all 
small e > 0. □ 

With this in mind, we subdivide integral (13.61) : suppose L roots meet on a set M. 
with codimA^ = £; without loss of generality, assume the coinciding roots are 
Ti(£), . . . , By continuity, there exists an e > such that only characteristic 

roots coinciding with r fc (£), k G {1, . . . , L}, in M. s are 7"i(£), . . . , Tl(^). Furthermore, 
we may assume that dM. £ G C 1 : for each e > there exists a set with C 1 boundary 
such that M. e c S e and meas(A1 £ ) — > meas(S' e ) as e — > 0. Then: 
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^1) Let Xm,e £ C°°(IR n ) be a smooth function identically 1 on A4 E and identically 
zero outside A4 2e ; now consider the subdivision of (13.61) : 



JB 2N (0) 

+ f e™<E^t){l-x M M)fXi)di; 
Jb 2N (o) 

for the second integral, simply repeat the above procedure around any root mul- 
tiplicities in B 2N (0) \M £ . 

(2) For the first integral, the case where the integrand is supported on Ai £ , split off 
the coinciding roots from the others: 

(3.7) f e^Eji^XMAOm)^ 

JB 2N (0) 

r L 

J B 2N (0) fc=l 

„ m 

+ / e <«*( Yl ^^(^XmAOHO^. 

JB 2N (0) k=L+l J 

(3) For the first integral, we use techniques discussed in Section [331 below to estimate 
it. 

(4) For the second there are two possibilities: firstly, two or more of the roots 
tl+i(0' • • • ; r m(0 coincide in M? e — in this case, repeat the procedure above for 
this integral. Alternatively, these roots are all distinct in Ai 2e — in this case, it 
suffices to study each integral separately as the Aj[,(£, t) are independent of t, and 
thus the expression (13.41) is valid and we can write 



„ m 

/ e<**( J2 e^A^t^XMAOHOd^ 

JB 2N (0) V ,„_r,i 7 

m „ 

£ / e^^A^OXMAOHOd^ 

._r , •, JBon(O) 



k=L+l 



k=L+l 



note that in this case we may use that the region is bounded to ensure the 
continuous functions are also bounded. 

Continue this procedure until all multiplicities are accounted for in this way. 

3.4.1. Roots separated from the real axis. The case where characteristic roots are 
separated from the real axis is similar to that for large |£|. Let us assume that Tfc(£) 
has no multiplicities in the set Q; now, a result similar to Proposition 13.31 holds for 
general integrals of this form, and thus 



Dr t D x(f e i(x<+Tm a(Of(0 dx 



< Ce- 5t 



LP 
Li 
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where p^ + q- 1 = 1, 1 < p < 2, N p > n(i - ±), r > 0, a a multi-index, / G C$°(R n ), 
5 > is a constant such that Imr(() > 5 for all £ G f2 and C = Cn,r,a, P > 0. So, in 
this case we have also have exponential decay of the solution. 

3.4.2. Roots meeting the real axis with finite order. In the case of bounded |£|, we 
must also consider the situation where the phase function r(£) meets the real axis. 
Suppose £o G is such a point, i.e. Imr(£o) — 0, while in each punctured ball 
around £o, B' e (£o) C Q, e > 0, there exists £ G -Bg(£o) so that Imr(£) > 0. Then, we 
claim that £o is a root of Imr(£) of finite order s: indeed, if £o were a zero of Imr(£) 
of infinite order, then, by the analyticity of Imr(£) at £o (which follows straight from 
the analyticity of r(£) at £o) it would be identically zero in a neighbourhood of £o> 
contradicting the assumption. 

In condition (H2) of Theorem 12.21 we actually have that there exist constants 
Co, c\ > such that, for all £ sufficiently close to £o, 

c \^-U s <\lmr(0\<Ci\^-U 2 . 
Indeed, the Taylor expansion of Imr(£) around £o, 

n 

lmr(0 = J>6 Imr(£ )te - (£„)i) + - £ | 2 ) , 
i=i 

is valid for £ G 5 e (£ ) C £1 for some small £ > 0. Now, if £ G B e (£ ), then -£ + 2£ G 
-B e (£o) also. However, 

n 

Imr(-£ + 2£ ) = - E ^ Imr (£o)(& - (&)*) + 0(|£ - £ | 2 ) ; 
i=i 

thus, for £ > chosen small enough, this means that either Imr(£) < or Imr(— £ + 
2£o) < — contradicting the hypothesis that Imr(£) > for all £ G f2; hence, 
^.Imr(( ) = for each z = l,...,n. In conclusion, Imr(£) = 0(|£ — £ | 2 ) f° r 
all £ G /4(£ ). 

Now, we need the following result, which is based in the calculation of the LP — L q 
decay estimate for the dissipative wave equation in [Mat 76] . but is here extended to 
a more general situation so that it can be used on a wider class of equations: 

Proposition 3.6. Let : U — > K, U C M™ open, be a continuous function and sup- 
pose £o G U such that 0(£o) — and that </>(£) > in a punctured open neighbourhood 
o/£o ; denoted by V \ {£o}- Furthermore, assume that, for some s > 0, there exists a 
constant Co > such that, for all £ G V, 

0(0 > co|£ -£o| s - 

Then, for any function a(£) that is bounded and compactly supported in U , and for 
all t > 0, / G L7 °°(M n ), and r G R, 

(3.8) / e -*0*|£ - &| r |a(0||/(OI ^ < C(l + , 

and 

(3-9) ||e-^-£oKO/(0|Uy) < C(l+t)-^||/|| x2 . 
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Proof. First, we give a straightforward result that is useful in proving each of the 
estimates: 

Lemma 3.7. For each p, M > and c > there exists C = C ps ^m,c > such that, 
for all t > 0, 

M 

x p e~ cxH dx < C(l + t)- {p+1)/ " and sup x p e~ cxH < C(l + t)- p/ " . 

0<x<M 

Proof. For < t < 1, each is clearly bounded: the first by M ^ and the second 
by M p . For £ > 1, set y — xt 1 ^; with this substitution, the first becomes 



while the second becomes 



y p r p "e- cy r 1 " dy < r (fi+1}/ " / y p e~ cy dy , 
Jo 



CO 



sup y p r p/ "e~ cyg < t~ p/ " sup ; 

(Ky^M* 1 /? 2/>0 

that the right-hand side of each is then bounded follows from standard results. □ 
Returning to the proof of (13.81) . as a(£) is bounded in U by assumption, we have 

/ c-*<«*ie-eorio(oii/(e)ide<c/ e-^ie-eoN/coide, 

</ v Jv 

where V = V H supp a; this, in turn, can be estimated in the following manner using 
the hypothesis on 0(£) and Holder's inequality: 

f e-^-U'\f(0\d£<C [ e-^-^-^\f(Z)\dt 
Jv Jv 

<C I e-^-^-Udafh-iV)- 
Jv 

Then, transforming to polar coordinates and using the Hausdorff- Young inequality, 
we find that, for some e > (chosen so that V C -E> e (£ ), possible since a(£) is 
compactly supported), 



e 



v 



dt\\f\\L~ { V) 



<C / l^-'e-^d^duWfW^), 



's™- 1 Jo 

Finally, by the first part of Lemma 13.71 we find 

-0(0* It _ £J r \n.(PM \ f(t\\rl£< a I n, r+n - l p- c oy st 



\t-U r Ht)\\m\di<C / f^-'e-^ dy\\f\\ Llm 
v Jo 

<C(l + ty^ s \\f\\ L i. 

This completes the proof of the first part. 
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Now let us look at the second part. By the second part of Lemma [3. 7\ 



< C{l + t) 



-2r/s 



-co\Z-Zo\ s t 



l/(0l 2 ^- 



The Holder inequality implies that 



e -coi^-eo[^|/ (e) | 2fi€ < sup | e - CQ |e-foi^|||/||2 <c\\f\\ 2 L2 , 

'V V 

and together these give the required estimate (I3.9p . 

So, using this proposition, we have, for all t > 0, and sufficiently small e > 0, 



□ 



did: 



a(o/(0<*e 



Be (Ho) 



£°°(R«) 



< / e 



-Imr(£)t 



o(0IK0l p ien/(0|d6<c(i + *)" B/ ' 



and, using the Plancherel Theorem, 



DfD° 



a(0/(0^ 



Be (&) 



L2(R S ) 



C||e <T( °V(O r ra(0/(0| 



^ 2 (^(«o)) 



<c||/IU" 



here we have used that |£l' a 'l r (0l r — ^ f° r ^6 ^' for r 6 N, a a multi-index. 
Thus, for all t > 0, 



3 i(:r-£+T(£)t), 



a(0/(0# , <C(l + f) 



:(!_!) 



'Be 

where 1 < p < 2, + = 1. 

Remark 3.8: If £o — 0, then Proposition 13.61 further tells us that 



D\D% 



a(0/(0# 



< 



L<J(R£) 



c(i + t) » 



LP 



If Rer(^o) = 0, then under condition (12.51) we have |r(£)| < |I mr (OI — c i|£ — £o| S1 
for £ near 6, and so we get 



did: 



a(0/(0# 



< 



L«(R») 



C(l + t)"V^Jl?-iJ 



L'' 



If both assumptions hold, we get the improvement from both cases, which is the 

_ I n+\a\+rs 1 \ I 1_ \\ 

estimate by C(l + t) \ s )\p i) . 
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3.5. Estimates for bounded |£| around multiplicities. Finally, let us turn to 
finding estimates for the first term of (13.71) . which we may write in the form 



r L 



where the characteristic roots 7i(£), . . . , tl(£) coincide on a set M. C Vt of codimen- 
sion £ (in the sense of Corollary 13. 5j) . O C M™ is a bounded open set and x £ C£°(Q). 

Unlike in the case away from multiplicities of characteristic roots, we have no ex- 
plicit representation for the coefficients (£, t), which in turn means we cannot split 
this into L separate integrals. To overcome this, we first show, in Section I3.5.1[ that 
a useful representation for the above integral exists that allows us to use techniques 
from earlier. Using this alternative representation, it is a simple matter to find esti- 
mates in the case where the image of the set M. is separated from the real axis. The 
argument may be extended to the case when it arises on the real axis as a result of 
all the roots meeting the axis with finite order. Such argument is more elaborate but 
not necessary for Theorem 12. 21 

3.5.1. Resolution of multiple roots. In this section, we find estimates for 

L 

]Te^4(£,t), 

k=l 

where Tl (£),... ,T£,(£) coincide on a set Ai of codimension £. For simplicity, first 
consider the simplest case, L = 2 and M. = {£o}; the general case works in a more 
involved but similar way. So, assume 

n(Co) = T2(6o) and T fc (£ ) ^ 7i(f ) for k = 3, . . . , m ; 

by continuity, there exists a ball of radius e > about 6o ; B £ (^o), in which the only 
root which coincides with Ti(£) is 72(6). Then: 

Lemma 3.9. For allt>0 and £ G B e (£ ), 
2 

(3.10) |^e^ (5) '^(£,t) < C(l + £) e -min(Imn«),ImT3«))* j 

fc=l 

where the minimum is taken over £ G B e (£o). 

Proof. First, note that in the set 

S := G R n : ^ r fc (0 VA; = 2, . . . ,m andr 2 (0 ^ r,(0 VZ = 3, . . . ,m} 

the formula (13 ,4p is valid for and ^(0- Now, recall that 

= Efc=i e iTfc «)*Aj (£, t) is the solution to the Cauchy problem (l3T2aD . (ErZcjl . 
and thus is continuous; therefore, for all 77 G W 1 such that 71(77) 7^ 7^(77) and 72(77) 7^ 
rfc(?7) for k — 3, . . . , m (but allow 7-1(77) = 7-2(17)), we have 

2 
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provided £ varies in the set S (thus, ensuring e in ®M.j(£) + e lT2 ^' t A^(^) is well- 
defined). Hence, to obtain (13.101) for all £ £ B e (£ ), it suffices to show 

|e <n ®*A-(£) + e i7 ^*A|(£)| < cte~ m[n( - lmTl( -^ ,lmT2( -^ t 

for all* > 0, £ £ B' e (£ ) = 5 e (£ ) \ {&}. 
Now, for all £ £ B£(£ ), * > °> 

(3.11) e iTl ^*Aj(0 + e iT2( ^^ 2 (0 

= sinh[(r 1 (0 - r 2 (£))t](e^^4(0 - ^^(0) 

+ cosh[(r 1 (0 - MZWK^AjiZ) + e^&Afc)) . 
Furthermore, we have the following estimates for all £ £ B' E (£ ), t > 0: 

(3.12) IsinhKnCO - r 2 (0)t](^(0e iT2( ^ - $(t)e**<&)\ 

< Ct\\e iT2 ^ t \ + |e in ^*|l < cte^ min( - lmT1 ^' liaT2 ^ t , 



(3.13) |cosh[(r 1 (0 - T2 {t))t){Ate)e'»® i + ^(fle 4 ^©*)] 

< Cte~ ia[n ( lmTl ^' liaT2 ^ t 

The proof of the first is simple: just note that 
sinh^O-r^))*] 



(Ti(fl- 73,(0) 



t as (ti(0-t 2 (O)^O 



or, equivalently, as £ — > £ through 5, and -<4j(£)(7i(£) — t 2 (£)) is continuous in 
.B £ (£o) for = 1, 2. The proof of the second is more technical and uses the explicit 
representation H3.4[) for the at points away from multiplicities of Tfc(£); otherwise 

it is similar and we omit it here. 

Combining (13.111) . (13 . 1 2[) and (13.131) we have (I3.10p . which completes the proof of 
the lemma. 

□ 

Suppose now that the characteristic roots tj. (£),..., Tl(£), 2 < L < m, coincide on 
a set M. of codimension £, and that 7~i(£) 7^ Tfe(0 for all £ £ M. when = L+l, . . . , m. 
By continuity, we may take e > so that the set .M £ = {£ £ M n : dist(£, M.) < e} 
contains no points 77 at which Ti(r]), . . . , Tl(t]) = 7^(77) for = L + 1, . . . , m. With 
this notation, we can extend Lemma 13.91 to the general situation: 



Lemma 3.10. For all t > and £ £ 7W e 

L 

|j^ e i7 *(f)*Aj(£,t) < C(l _|_t)i-l e -*mm fc=1 ,..., L Im^(S) 
fc=l 

where the minimum is taken over £ £ B E (^ Q ). 



Note that this estimate does not depend on the codimension of A4 nor its geometric 
structure. 
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3.5.2. Phase function separated from the real axis. We now turn back to finding LP — 
L q estimates for 

r L 

J Q elx< ( E ^ (0t 4(£> *)) xiom # , 

when Ti (£),... ,tl(£) coincide on a set M. of codimension £; choose e > so that 
these roots do not intersect with any of the roots tl+i(£), . . . , T m (£) in .M e . 

Here we can assume that there exists 5 > such that Imr^O > 5 for all £ G .M £ — 
so, min fc IniTfc(£) > 5. For this, we use the same approach as in Section [3.4.11 but 
using Lemma [3. 101 to estimate the sum. Firstly, the L 1 — L°° estimate: 



f L 



n k=x 



<maxsup|r fc (0| r / YV^U^t) |f |H|/(0| <fc 



n 1 



< c(i + t) L - 1 e - <5< 

Similarly, the L? — L? estimate: 



k=X 

L°°(n)<C(l+t) L - L e 



L-l-St 



C L 

D r t D a x ( je*< ( e^A)^ f ))x(fl/(fl 



fe=i 



^ e ^^4(e,t)r fc (o r )rx(o/(0 



fe=l 



£-1 -5t 



< C(l + t) L - L e 



L 2 {Q) 



L 2 (n) 

L-X-St 



< C(l+t) L - L e 



L 2 



Then, 
Dl 



r L 

D * ( J n ^ ( E *))x(0/(£) dx 



L-X-St 



<C(l + t)- ~e 

where p~ l + q~ x = 1, 1 < p < 2. Once again, we have exponential decay. 



LP 



4. Applications 

In this section we will briefly consider applications of Theorems 12.11 and 12.21 to 
Fokker-Planck equations and wave equations with dissipation and negative mass. 
There are further applications to Grad systems linearised near equilibrium points 
where Theorems 12.21 and 14.11 immediately yield the corresponding decay rates. The 
size of these systems depends on the number of moments and dimension of the space. 
Some examples of these systems and their stability has been analysed in [VR03| . 
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4.1. Fokker— Planck Equation. The classical Boltzmann equation for the particle 
distribution function / = /(£, x, c), where i,c6 W 1 , n = 1,2, 3, is 

(d t + c-V x )f = S(f), 

where S(f) is the so-called integral of collisions. The important special case of this 
equation is the Fokker-Planck equation for the distribution of Brownian particles, 
when the integral of collisions is linear and is given by 

n 

S(f) = V c ■ (C + V c )/ = d cM + d Ck )f. 

k=l 

In this case the kinetic Fokker-Planck equations takes the form 

f n \ n 

dt + Yl ° kd ^ ) x ' c ) = Yl dc ^ Ck + dc *)f- 

k=l J k=l 

The Hermite-Grad method of dealing with Fokker-Planck equation consists in de- 
composing f(t,x, •) in the Hermite basis, i.e. writing 



f{t,x,c) = —m a (x,t)ip c 



|«|>o 

where ip a (c) = (27r)~ n / 2 (— <9 c ) a exp(— ^j-) are Hermite functions. They are derivatives 
of the Maxwell distribution ip° which annihilates the integral of collisions and form 
a complete orthonormal basis in the weighted Hilbert space L^(R n ) with weight 
w = l/ip°. This decomposition yields the infinite system 

dtmp(x,t) + (3 k d Xk m/3- ek (x,t) + d Xk mi3 +ek (x,t) + \P\mp(x,t) = 0. 

The Galerkin approximation / of the solution / is 



f N (t,x,c) = —^m a (x,t)ij c 



0<\a\<N 



with m(x,t) = {mp(x,t) : < |/3| < N} being the unknown function of coefficients. 
For m(x, t) one obtains the following system of equations 

D t m(x, t) + AjD Xj m(x, t) — iBm(x, t) = 0, 
i 

where B is a diagonal matrix, B a $ = \a\8 ai p, and the only non-zero elements of the 
matrix Aj are aj e -"° = ay, a" +ej,a = 1. For details of these calculations see [VR04] . 
Hence, the dispersion equation for the system is 

P(r, = det(r/ + £ - iB) = 

(4.1) " N 

P(r, 0) = det(r/ - iB) = r JJ(r - ji)^ . 

Properties of this polynomial P(r, £) have been extensively studied by Volevich and 
Radkevich in [VR04] . who gave conditions and examples of situations when Imr,-^) > 
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0, for all £ 7^ 0. In our situation here we have to take additional care of possible 
multiple roots, as is done in Theorem 12.21 

Assume now that P(t, £) is a stable polynomial, i.e. its roots r(£) satisfy Imr(£) > 
and Imr(£) = imply £ = 0. We will say that P(r, £) is strongly stable if, 
moreover, its roots r(£) satisfy Imr(() > e > for large £. It follows that we satisfy 
the conditions of Theorem 12.21 and we have to determine the order with which the 
characteristic arrives at the origin. We have the following theorem about time decay 
of solutions to Cauchy problems for equations with strongly stable symbols. 

Theorem 4.1. Let a strongly stable polynomial P(r, £) of order m have a strictly 
hyperbolic principal part and assume that d T P(0,0) 7^ 0. Let a be the multiindex of 
the smallest length \a\ such that <9|P(0,0) 7^ 0. Let u(x,t) be the solution of the 
Cauchy problem P(D t , D x )u = 0,d l t u\ t =o = fi, < I < m — 1. Let 1 < p < 2, 
2 < q < 00, and 1/p + 1/q = 1. Then 



m—l 



\ U (;t)\\L,<C(l+t)-m- p - q )Y^\\fl\\ wNl , 



1=0 



where iVj = ^— £([n/2] + 1) — I. Moreover, we have the estimate 



m—l 



dld?M;t) <C(l+t)-^ +r )(^)j2Wf' 

v x) 1=0 



Indeed, since the polynomial is strongly stable, the estimate in Theorem 14.11 follows 
from Theorem 12.21 and Remark 12.31 The improvement in the last estimate for deriva- 
tives comes from the fact that there is only one root r such that r(0) = and so the 
last statement of Remark 12.31 applies. In certain cases it can be shown that actually 
I ct I = 2, in which case we have the same decay as for dissipative wave equation. 

Let us write the polynomial P(r, £) from (14.1 j) in the form 

M+l 
j=0 

where Pj is a homogeneous polynomial of order M+l — j, and assume that P(t, £) is 
strongly stable. In |VR04j it was shown that Pm+i(t,£) = 0, 7mPm(TjO = M\t for 
some 7 M > 0, and 7m-iPm-i(t,0 = Ml ih^ ~ M ^\ 2 for some 7m-i > 0. 

It can be readily verified now that conditions of Theorem 14 . 1 1 hold with \a\ = 2, from 

which we get the estimate with (1 + t) 2[ -p *> . 



4.2. Wave type equations with (negative) mass and dissipation. Here we 
will show that we can still have time decay of solutions if we allow the negative mass 
but exclude certain low frequencies for Cauchy data. This is given in (14. 2[) below. 
Nonnegative but time dependent mass and dissipation with oscillations have been 
considered before. See, for example, 11 1>()3] and references therein. 
Let us consider second order equations of the following form 

d\u — c 2 Au + 5d t u + fiu = , 
u(0,x) = 0, ut(0,x) = g(x) . 
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Here 8 is the dissipation and \i is the mass. For simplicity, the first Cauchy data 
is taken to be zero. The general case can be treated in the same way. Let us now 
apply Theorem 12.11 to the analysis of this equation. The associated characteristic 
polynomial is 

T 2 -c 2 \^-i8r-n = 0, 

which has roots 

r±(0 = y ± v/ c 2|e + /i-574. 
Now, we have the following cases, which correspond to different cases of Theorem 

EU 

• 8 = fi = 0. This is the wave equation. 

• 8 — 0, [A > 0. This is the Klein-Gordon equation. 

• /J, — 0, 8 > 0. This is the dissipative wave equation. 

• 8 < 0. In this case, Imr_(£) < | < for all £, hence we cannot expect any decay 
in general. 

• 8 > 0, /i > 0. In this case the discriminant is always strictly greater than — 8 2 /4, 
and thus the roots always lie in the upper half plane and are separated from the 
real axis. So we have exponential decay. 

Here is the main case for us: 

• 8 > 0, ji < 0. In this case, note that Imr_(£) > if and only if c 2 |£| 2 + \i > 0, i.e. 
the critical value is |£| = ^/\JT\/c. Therefore, the answer depends on the Cauchy 
data g. In particular, if supp g is contained in c |£| 2 + /i > 0, then we may get decay 
of some type. More precisely, let -8(0, r) be the open ball with radius r centred at 
the origin. Then we have: 



if g is such that supp PI -6(0, ^-M) -l: ; then we have no decay; 



if there is some e > such that supp g C lR n \ -8(0, -~r- then the roots are 
either separated from the real axis (if 8 > 0), and we get exponential decay, or 
lie on the real axis (if 5 = 0), and we get Klein-Gordon type behaviour (since 
the Hessian of r is nonsingular). 



if, for all g, supp g C R n \ B(0, ^) = <^ |f | > ^ } , then again we must 



consider 8 = and 8 > separately. 

If 8 = 0, then the roots lie completely on the real axis, and they meet on the 
sphere |£| = ^/J/xj/c. It follows from (12.41) with L = 2 and t = 1 that, although 
the representation of solution as a sum of Fourier integrals breaks down at the 
sphere, the solution is still bounded in a (1 /^-neighbourhood of the sphere. 
In its complement we can get the decay. 



If 8 > 0, then the root r_ comes to the real axis at |£| = j n which case 

we get the decay 

(4.2) IM,*)l|w < C{l + t)-(p-^\\g\\ LP . 

Indeed, in this case the order of the root r_ at the axis is one, i.e. estimate 
(TOD holds with 3 = 1. Here 1/p + 1/q = 1 and 1 < p < 2. Note also 
that compared to the case of no mass when £ = n, the codimension of the 
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8et ^ e R» : |<| = XM} is I = 1. We can apply the .a 8 t ease of Part II of 
Theorem 12.11 with L = 1 and s = £ = 1 which gives estimate (14. 2 j) . 
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